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Abstract 

Procedures for conecting the beam-beam effects in luminosity measurement at 
CLIC at 3 TeV CM energy are described and tested using Monte Carlo simulations: 

• Correction of the angular counting loss due to the combined Beamstrahlung 
and initial-state radiation (ISR) effects, based on the reconstructed velocity of 
the collision frame of the Bhabha scattering. 

• Deconvolution of the luminosity spectrum distortion due to the ISR emission. 

• Correction of the counting bias due to the finite calorimeter energy resolution. 

All procedures were tested by simulation. Bhabha events were generated using 
BHLUMI, and used in Guinea-PIG to simulate the outgoing momenta of Bhabha 
particles in the bunch collisions at CLIC. Residual uncertainties after correction are 
listed in a table in the conclusions. The beam-beam related systematic counting 
uncertainty in the luminosity peak can be reduced to the order of permille. 
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23 1 Introduction 



24 Luminosity, L, and luminosity spectrum, ^{Ecm), are key input to analyses of most mea- 

25 surements at a linear- collider, including mass and cross-section measurements, as well as the 

26 production-threshold scans. Precision measurement of the luminosity is thus essential for the 

27 physics programme at a linear collider. The standard way to measure luminosity is to count 

28 Bhabha-scattering events recognized by coincident detection of showers in the fiducial volume 

29 (FV) in both halves of the luminometer in the very foward region in a given energy range. The 

30 number of events A'^ is then divided by the Bhabha cross section a integrated in the corresponding 

31 region of the phase space. This can be formally expressed as follows, 

a(Z(aff,£f;f))' 

32 Here E{Q.'"2,E["2) is a function describing the selection criteria for counting the detected 

33 events, and Z(Q^^, Sp^) is a function describing the corresponding region of phase space where 

34 the cross section is integrated. These functions can be expressed as products E = Ylj and Z = 

35 rTi Q where the functions and Q are based on specific topological and kinematical properties 

36 of the detected/generated pair. For each /, the physical meaning of and Q con^esponds to 

37 each other, although their mathematical fonii may be different The set of functions and 

38 i^i includes the angular selection requiring both particles to be detected in the FV, as well as the 

39 energy range selection and possible further cuts to eliminate background. 

40 Bhabha scattering is a well known QED process, for which cross-section calculations with 

41 relative uncertainty better than 10^^ are available [1]. However, at high beam power, the energies 

42 and the polar angles of the Bhabha particles are additionally strongly influenced by the beam- 

43 beam effects [2, 3]. Because of the random and asymmetric momentum loss when electrons 

44 emit Beamstrahlung, the CM frame of the Bhabha process moves with respect to the lab frame 

45 with axial velocity different for every colliding pair. As a consequence, S and Z operate on 

46 kinematical arguments in different reference frames. Thus, if E and Z have the same form, 

47 different regions of the phase space will be covered, leading to a systematic bias in the luminosity 

48 measurement. This systematic bias cannot be neglected at the future linear colliders due to the 

49 high energy and charge density, and is particularly accute at the 3 TeV CLIC. 

50 A way around this problem is to define S and Z such that the counting rate is independent of 

51 the reference frame. Some of the functions and Q can be defined invariant to the boost along 

52 the beam axis. This is, for example, the case with the cuts on the reconstructed CM energy. 

53 However, the requirement that the outgoing particles hit the FV of the detector on both sides 

54 does not posess such invariance. In this paper, a definition of ^py and i^fv is proposed such that 

55 both the experimental count and the cross-section a are reconstructed in the same reference 

56 frame, namely the collision frame, which will be properly defined in the appropriate section. 

57 The physical processes affecting the luminosity measurement will be outlined and the used 

58 terms and notation defined in Sec. 2. The analysis method with the correction procedures, as 

59 well as the test results will be described in Sec. 3. In the conclusions, the main advantages of 

60 the presented method are restated, and the final uncertainties are listed and briefly discussed. 



See in particular Sec. 3.2 and Eqs. 6 and 7 
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61 2 The physical processes affecting the luminosity 

62 measurement and the outline of the correction procedure 

63 The sequence of physical processes relevant to the present discussion is schematically repre- 

64 sented in Fig. 1. Due to the pinch effect during the bunch collision, both particles may emit 

65 Beamstrahlung photons and so lose energy and momentum before the interaction. Thus in gen- 

66 eral, Ecm < Eq = lEheam- The CM energy distribution at this stage is the actual luminosity 

67 spectrum ^{Ecm)- The probability of the Bhabha scattering scales with 1/^=1 /-Ecw resulting 

68 in the CM energy distribution of the Bhabha events .'^{Ecm) °^ -^{Ecm) /^cm- Bhabha pro- 

69 cess is itself accompanied by emission of the initial-state radiation (ISR) that is nearly collinear 

70 with the initial particle momenta, as well as the final-state radiation (FSR) that is approximately 

71 collinear with the outgoing particle momenta. Since the ISR is nearly collineai^ with the beam 

72 axis, it misses the luminometer, so that the CM energy reconstructed from the detected particles 

73 is EcM; rec < ^CMf 'ind the Corresponding spectrum is, 

^max 

KEcM.rec)= I I^{Ecm)-^^ {^^^)dEcM (2) 

J Ecm Ecm 



74 where J^{x) is the distribution of the fractional CM energy losses due to ISR. ^{x) is approx- 

75 imately independent of Ecm- 

76 Due to the finite energy resolution of the LumiCal, the reconstructed spectrum is smeared, 

77 which can be represented as a convolution with a normalized Gaussian-' . 




Beamstrahlung ISR FSR Synchrotron 

radiation 



Figure 1: Schematic representation of the physical processes affecting the luminosity measure- 
ment 



^Strictly speaking, the smearing width depends on the deposited energy of the showers. However, as only a relatively 
narrow energy range is being analyzed here, the smearing width will be treated as being approximately constant. 
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78 The terni collision frame will be used here for the frame of the two-electron system^' after 

79 emission of Beamstrahlung and ISR and before emission of the FSR"*'. The scattering angle in 

80 the collision frame is denoted B'^"'' . Due to the radiation prior to the colUsion, the collision frame 

81 has a non-zero velocity ^coU, and the outgoing particle angles in the lab frame, and 62'^, 

82 are not symmetric. In a significant fraction of the total number of events, the acoUinearity is so 

83 large that the two particles aix not detected in coincidence within the FV of the LumiCal. In this 

84 way, Beamstrahlung and ISR induce an angular counting loss of Bhabha events. 

85 Finally the electromagnetic deflection (EMD) of the outgoing electrons in the field of the 

86 opposing bunch induces a small additional angular counting loss. 

87 The outline of the procedure of the Bhabha-count analysis is as follows: 

88 1. Reconstmct the CM energy Ecm.rec and the collision-frame velocity ^coii for each pair 

89 detected in the FV of the LumiCal, from the angles and the measured particle energies. 

90 2. Assign weight to events to correct for reduced acceptance due to ^coiu as shown in Sec. 

91 3.2. 

92 3. Deconvolution of the ISR energy loss ^{x) from the spectrum h* {EcM.rec), in order to 

93 restore the 3S*{Ecm) CM energy spectrum of the Bhabha events (see Sec. 3.3). 

94 4. Integrate 3S*{Ecm) in the energy range of measurement. 

95 5. Correct the systematic effect of the finite energy resolution of the LumiCal on the number 

96 of counts in the peak (Sec. 3.4). 

97 The absolute luminosity in the measured energy range is then given by Eq. 1 , and the approx- 

98 imate differential form of the luminosity spectmm with the LumiCal energy smearing can be 

99 obtained as ^* (Ecm) = SB* {Ecm)EIm- 



100 3 Analysis and correction procedures 

101 3.1 Simulation methods used to test the analysis procedure 

102 To test the analysis procedure, Bhabha events in the bunch-collision were simulated with the 

103 Guinea-PIG software [3]. The initial bunch coordinate- and momentum distributions were taken 

104 from the simulation results by D. Schulte et al. [4]. The coordinate distribution covered more 

105 than 10 a bunch widths both in the horizontal and the vertical directions. The angular distribu- 

106 tion of the particles in the bunch was quasi-Gaussian, with transverse emittance of 660 nm rad in 

107 the horizontal, and 20 nm rad in the vertical direction. The bunch collision was simulated in the 

108 CM frame of the colliding bunches, which is equivalent to a head-on collision with zero cross- 

109 ing angle. The beam overlap reduction due to the crossing angle is offset by the crab-crossing 

110 scheme. 

■^Unless stated otherwise, electron always refers to electron or positron 

^In reality, ISR and FSR can not be cleanly separated even theoretically, due to the quantum interference between 
them. Thus in practice, the collision frame is defined as the CM frame of the final electrons together with all 
radiation within a given tollerance angle with respect to the final electron momenta. The assumption of clean 
separation between ISR and FSR introduces a small uncertainty in the final result. 
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111 The Bhabha events were produced using a method resembling that by C. Rimbault et al. [5]: 

112 After generating the initial four-momenta of the colliding pairs, the decision is made by 

113 Guinea-PIG whether the Bhabha scattering will be realized in the collision, based on the \/s 

114 proportionality of the Bhabha cross section. If a Bhabha event is to be realized, the final four- 

115 momenta are picked from a file generated at 3 TeV by the BHLUMI generator [6]. The momenta 

116 are then scaled to the CM energy of the colliding pair, rotated to match the collision axis, and 

117 then boosted back to the lab frame. Finally the outgoing Bhabha electrons are tracked to simulate 

118 the electromagnetic deflection. Nearly four million Bhabha events were generated in the polar 

119 angle region between 10 and 200 mrad in the generator frame, with additional cuts between 37 

120 and 90 miad in the collision frame. 

121 The interaction with the detector was approximated in the following way: 

122 • The four-momenta of all electrons and photons within 5 mrad of the most energetic shower 

123 were summed together on each side. The 5 mrad criterion corresponds closely to the 

124 Moliere radius of the high-energy showers in the LumiCal [7]. The Beamstrahlung pho- 

125 tons were not included as they are emitted close to the beam axis. For synchrotron radi- 

126 ation, the chai^acteristic emission angles are of the order I/7, which is smaller than 10^^ 

127 mrad for electron energies in the TeV range, therefore photon and electron four- vectors 

128 can be added. 

129 • The energy resolution of the LumiCal was included by adding random fluctuations to 

130 the final particle energies. The random fluctuations were sampled from the Gaussian 

131 distribution with an energy-dependent standard deviation Oe = a\fE, with a = 0.21 [8, 9]. 

132 • The finite angular resolution of the LumiCal was included by adding random fluctuations 

133 to the final particle polar angles. The nominal value of Oq = 2.2 x 10^^ estimated for the 

134 ILC version of LumiCal [7] was used. Higher values for Oq were also tested. 

135 The FV of the LumiCal at CLIC corresponds to the angular range between 43 and 80 mrad 

136 around the respective beam axes. 



137 



3.2 Invariant counting in the collision frame 



138 The movement of the collision frame with respect to the lab frame is responsible for the acollineai- 

139 ity leading to the angular counting loss. The velocity of the collision frame with respect to the 

140 lab frame j3co//, can be approximately calculated from the measured polar angles. If ^coU is taken 

141 to be coUinear with the z-axis, the expressions for the boost of the Bhabha scattering angles into 

142 the lab frame give, 

_ sin(er + e^-^) 

'^'^"""sinef^^ + sinef^ ^ ^ 

143 Eq. 4 does not depend on any assumptions about the number of emitted ISR and Beam- 

144 strahlung photons, nor on their direction, apart from the assumption that the vector sum of their 
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145 momenta is collinear with the z-axis^^ . 

146 If events from a subset characterized by a given Pcoii ai"e plotted in the | tan 02 1 vs. |tan6i| 

147 graph, they he on a curve displaced from the central diagonal, as schematically represented by 

148 the dashed hne in Fig. 2. As can be seen from the figure, the range of accepted scattering angles 

149 decreases with increasing Pcoii- The effective limiting angles 6^^^^^ and d^^'H for the subset of 

150 events charaterized by a given PcoU ai"e obtained by boosting ^i^n and Omax into the collision 

151 frame. 




Figure 2: Schematic representation of the distortion of the polar angles due to the movement of 
the collision frame. The box represents the region in which both electrons hit the FV, 
and the dashed line represents the event subset characterized by a given jicoii ■ Q^^'l and 
^max denote the effective limiting scattering angles for this subset. 



152 To account for the smaller acceptance of the events characterized by a given ^coiu every event 

153 has to be weighted with the appropriate correction factor. In this way, the number of events 

154 between ^min and O^ax in the collision frame is recovered for each ^coU subset separately. The 

155 weighting factor is defined in the following way: 



^ Strictly speaking, ^coU has a small radial component j8p, which is larger than 0.01 in only about 5 permille of 
cases. However, the influence of /3p on the polar angles of the Bhabha pair is almost indistinguishable from 
an additional axial boost. Thus for the purpose of recovering the counting loss due to the acollinearity, ficoll is 
approximately treated as a scalar quantity. 
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^iPcoi,) = '-^, . (5) 

ecoll 
min 



156 The FV selection function is thus defined as^^ 

w ; 

; otherwise 



:,FV 



157 Using this FV selection function, the number of events N satisfying the condition 0"'" G 

158 (0mini ^max) in the Collision frame is reconstructed. The coiTcsponding function i^py for the 

159 cross-section integration is thus, 

_ f 1 ; S'^"" £ (^min, ^max) 

1 ;otherwise 



3.2.1 Test of the collision-frame counting method 



161 To test the counting method, histograms of Ecm.rec reconstructed from kinematic parameters of 

162 the detected particles were generated in the following way: 

163 Control histogram : All events with the scattering angle in the collision frame 0"'" such that 

^min < < ^max are accepted. Therefore this histogram is not affected by counting 

165 losses due to Beamstrahlung and ISR. This is, of course, only possible in the simulation. 

166 Uncorrected histogram : Events hitting the FV of the LumiCal in the lab frame. 

167 Corrected histogram : Events hitting the FV of the LumiCal in the lab frame, stored with 

168 the weight w calculated according to Eq. 5 

169 The full kinematical information, including the particle energy, was used for the reconstruc- 

170 tion of the CM energy. The four-momenta of coUineai^ photons were added to the electron 

171 four-momenta as described in Sec. 3.1. Apart from the coUinearity criterion, no distinction was 

172 made between the ISR and the FSR photons. 

173 To calculate the correction weight w, the approximate expression for the angular differential 

174 cross section do /dd ^6^^ was used. 

175 The results are shown in Fig. 3. The control spectrum is plotted in black, red is the spectrum 

176 affected by the counting loss, green is the coixected spectrum. 

177 The blue line represents the events inaccessible to the correction due to the high value of ^coii- 

178 In the subsets of events characterized by ^coii above a certain treshold, at least one electron is 

179 always lost (see Fig. 2). However, for such events, the Beamstrahlung-ISR energy loss is also 



^By the standard definition of the polar angle 6, the interval corresponding to the FV on side of the IP is (6mini ^max) 
and on the other side of the IP, [tz — Omnx > ^ — dmin) 
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180 above a certain minimum, so that they are only present in significant number below 2200 GeV. 

181 A small number of high-/3c„// events ai^e also present at energies above 2200 GeV, as seen in the 

182 zoomed figure (Fig. 3, right), where these events are scaled by a factor 100. In these events, 

183 pcoH has a relatively high radial component, due to the off-axis radiation before colhsion. This 

184 increases the acollinearity of such events relative to other events with similar energy loss (see 

185 footnote 5). The relative contribution of these events to the peak integral above 95% of the 

186 nominal CM energy is of the order of 2 x 10^^. 

187 Before correction, the counting loss in the peak integral above 95% of the nominal CM energy 

188 was 3.8%. After correction, the residual relative deviation in the peak with respect to the control 

189 spectrum is (-0.1 ±0.4(stat.)) x 10"^. In the tail between 80% and 90% of the nominal CM 
energy, the counting loss before connection was 43.1%. After coiTection, the residual relative 



190 



191 deviation in the tail is (—3.6 it 1.8(stat.)) x 10^ , which includes a deviation of (—2.7 ±0.1) x 

192 10^^ due to the lost events in the tail. The statistical uncertainty of the residual deviation was 

193 estimated taking into account the correlations between the con^ected and the control spectra. The 

194 precision of the Beamstrahlung-ISR correction is of the order of permille despite the presence 

195 of the following sources of systematic uncertainty of the correction: 

196 • The assumption that the deformation of the Bhabha angles induced by Beamstrahlung and 

197 ISR is well described as a Lorentz boost along the beam axis (this is the source of the 

198 "lost" events in the peak), 

199 • The use of the approximate angular differential cross section for the Bhabha scattering in 

200 the calculation of w, 

201 • Assumption that all ISR is lost, and all FSR is detected (this assumption has, in principle, 

202 an influence on the calculation of pcoii, and consequently of w). 




EcM (GeV) (GeV) 



Figure 3: Correction of the counting loss due to Beamstrahlung and ISR. Left: whole spectrum; 

right: zoom on energies above 2200 GeV. Black: Simulated control spectrum with- 
out counting loss due to Beamstrahlung and ISR; red: Reconstructed Ecm spectrum 
affected by the counting loss; green: Reconstructed spectrum with correction for the 
counting loss due to Beamstrahlung and ISR; blue: events inaccessible to the correc- 
tion due to high ficoii (see text). 
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Table 1 : Dependence of the uncorrected and the corrected peak-count deviations on the polar- 
angle uncertainty of the LumiCal 



Ad (mrad) Uncoixected (10^^) Corrected (10^^) 



0.010 -38.1±0.3 -0.0±0.4 

0.022 -38.2±0.3 -0.1±0.4 

0.05 -38.3±0.3 -0.1±0.4 

0.10 -38.8±0.3 -0.3±0.4 

0.20 -40.1±0.4 -0.4±0.4 

0.50 -46.8±0.4 -1.6±0.4 



203 A test of this method for the case of ILC at nominal CM energies of 0.5 and 1 TeV, and for a 

204 wide range of beam-paiameter uncertainties is reported in Ref. [10]. After coiTection of the lost 

205 events, precision of a fraction of permille is achieved in the upper 80% of energy, independently 

206 of the precision with which the beam parameters are known. 

207 3.2.2 Effect of the particle energy- and polar-angle resolution on the 

208 Beamstrahlung-ISR counting loss 

209 The relative energy resolution of a calorimeter can be parametrized as a quadratic sum of an 

210 energy-dependent stochastic term and a constant shower leakage term [11]. 



21 1 The 40-layer design of the LumiCal for CLIC is assumed to be well described by the paramet- 

212 rization with the constant term b equal to zero [8]. A conservative limit for the constant term can 

213 be taken to be 1.1%, based on the analysis in Ref. [11], made for the ILC version of LumiCal, 

214 with only 30 layers (implying considerable leakage at 1.5 TeV). The correction procedure was 

215 therefore tested with three different values of the constant term b: 0, 0.35% and 1.1%, while the 

216 stochastic term was kept at a = 0.21. In all these tests, both the uncorrected Beamstrahlung-ISR 

217 counting loss, and the residual deviation after coiTection agree within their respective statistical 

218 uncertainties. Only results for Z? = are presented in this paper 

219 Tests were also performed with different values of the polar-angle uncertainty of the LumiCal 

220 in the range between 0.01 mrad and 0.5 mrad, beside the nominal value of 0.022 mrad [7]. The 

221 results of these tests are shown in Tab. 1. One may note that the uncoixected counting loss 

222 increases significantly when the polar-angle uncertainty increases above 0. 1 mrad. The effect on 

223 the count deviation in the corrected peak is, however, about six times smaller 

224 3.3 Deconvolution of the ISR energy loss 

225 After correcting for the angular counting loss, the ISR energy loss can be deconvoluted from 

226 the resulting spectrum h{EcM.rec) to restore the CM energy spectrum of the Bhabha events 




(8) 
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227 ^*{EcMy\ When Eq. 2 is binned 

228 discrete fomi, 



in N sufficiendy naixow bins, it takes approximately the 



N 



(9) 



229 



Where, 



hi = 



iAE 

/ h{EcM,rec)dEcM,. 



,rec 



ii-l)AE 
JAE 




'•-1/2 
j-i 



(10) 



i-l/2 



230 As the matrix has triangular form, eq. 9 can be solved for exactly, using the Jacobi 

231 method. The solution proceeds from high-energy towards the lower-energy bins, indroducing an 

232 increasing uncertainty towards lower energies. 

233 To obtain J^^-, the function was parametrized by the beta distribution used for the 

234 parametrization of the beam spectra of linear colliders [12], 



235 The parameters were obtained by fitting Eq. 1 1 to the fractional CM energy distribution 

236 after ISR, reconstructed from the same BHLUMI data set as used in Guinea-PIG. The fit was 

237 performed with variable binning in order to have sufficiently fine binning near x = I, while 

238 avoiding large differences in statistical uncertainties for individual bins. The data histogram 

239 was first normalized to the unit integral. The results are shown in Fig. 4. The parameter 

240 was obtained as the ratio of the number of counts in the narrow peak above x = 0.99995 to the 

241 number of counts in the entire spectrum, and the remaining coefficients were obtained by fitting 

242 the function to the data in the range (0.7, 0.99995) ^) 

'The Bhabha-event spectrum is marked with a star here, because it is smeared by the finite energy resolution of the 
LumiCal. See Sec. 3.4. 

^ The angular cuts in the lab frame cause significant losses in the distribution for x < 0.5 because high energy loss in 
ISR emission correlates with high acollinearity. This affects the overall normalization, and thus the value of oq. 
The value of obtained here is appropriate for the deconvolution of the simulated spectrum where the same set 
of BHWIDE samples was used. However, for the analysis of the real experimental data, the distribution without 
cuts in the lab frame should be used. 

''The functional form of Eq. 1 1 suggests that the ratio ao/ai can be fixed by the normalization requirement. How- 
ever, the beta distribution fails to properly describe the form of J^(x) for x < 0.7 (regardless of the angular cuts 




(11) 
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Fractional ISR energy loss 
30 = 0.31677 ±0.00032 
8^ = 0.1331 ±0.00049 

82 = -0.364 ±0.024 

83 = -0.89374 ±0.00056 



0.8 0.85 0.9 0.95 1 

X 

Figure 4: Fit of the relative energy-loss distribution due to the ISR. 

243 From the normalization of the data histogram, it follows that, ao + Cai = 1, from which 

244 the covariance of ao and ai was calculated as cov(<3o,fl'i) = {{ao— < oq >){ai — < ai >)) = 

245 — (1/C) {{ao— < ao >)^) = — (l/C)(7o = —ai/{l —ao)OQ, where ao is the statistical uncerat- 

246 inty of the parameter ao- The correlation of ao with ^2 and 03 was assumed to be zero, and the 

247 correlation coefficients between the fitted parameters ai 2,3 were calculated from the covariance 

248 matrix obtained in the fit procedure. The nonzero correlation coefficients are the following: 

249 po,i = -13%, pi,2 = 83%, pi,3 = 96%, P2,3 = 75%. 

250 3.3.1 Test of the deconvolution procedure 

251 In this test, the following histograms were generated: 

252 Control histogram was filled with simulated CM energies before ISR emission, and then 

253 smeared with a normalized Gaussian with constant width corresponding to the LumiCal 

254 energy-resolution at the peak energy. 

255 Histogram with ISR energy loss h{EcM,rec) is the same as the control histogram from Sec. 

256 3.2 - filled with energies reconstructed from the final-state kinematics, and with inclusion 

257 of the LumiCal energy resolution. 

258 Deconvoluted histogram was obtained by solving the system of linear equations repre- 

259 sented by Eq. 9, taking the binned data of the affected histogram as hj. 

260 For each histogram, event selection was made on the scattering angles in the collision frame, 

261 so that the Beamstrahlung-ISR angular counting loss is not present. This was done in order to 

262 assess the accuracy of the deconvolution separately from the Beamstrahlung-ISR counting-loss 

263 correction. Results are shown in Fig. 5. 

264 Before deconvolution, the relative counting loss in the peak above 95% of the nominal CM 

265 energy was 23.4%. After deconvolution, the relative residual deviation of the peak integral with 

in the lab frame discussed above), so that the overall norm is different than the integral of the beta distribution 
extrapolated from the fit. Therefore, a\ was allowed to vary freely in the fit. 



0IO 



lO'' 
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2200 2400 2600 2800 3000 

EcM (GeV) 

Figure 5: Deconvolution of the ISR deformation of the luminosity spectrum. Yellow: the con- 
trol histogram - simulated Ecm before emission of ISR, smeared with a normalized 
Gaussian; black: the histogram affected by the ISR energy loss - reconstructed Ecm 
from the detected showers, green: deconvoluted spectrum 

266 respect to the control histograms is (+1.3 ±2.1) x 10^^. In the tail between 80% and 90% of 

267 the nominal CM energy, the ISR energy loss increases the count by 14.5%. After deconvolution, 

268 the remaining deviation in the tail is (—2.3 it 3.9) x 10^^. 

269 The contributions from the uncertainties of the fitted parameters of the ISR energy-loss func- 

270 tion '^{x) were added to the statistical uncertainty of the residual deviation after deconvolution. 

271 The full covariance matrix of the fit parameters was used, together with the partial derivatives of 

272 the count estimated by variation of the fit parameters by one sigma, one parameter at a time. With 

273 the statistic of about four million generated Bhabha events, the uncertainties due to the fit pa- 

274 rameters are {5N /N)peakjSRfit = 0.53 x 10"^ for the peak, and {5N /N)taUjSRfit = 0.07 x 10"^. 

275 When the deconvolution step is tested with a non-zero leakage term in the LumiCal energy 

276 resolution, the residual uncertainties agree with the ones presented above within the statistical 

277 uncertainties. 

278 3.4 Effect of the LumiCal energy resolution on the counting rate in the 

279 peak 

280 The finite energy resolution of the LumiCal introduces a counting bias in two ways: 

281 1. By asymmetric redistribution of events from each side of the sharp energy cut E^ut used to 

282 define the peak, due to the slope of the underlying distribution at the position of the cut. 

283 2. By cutting off a portion of the low-energy tail of the quasi-Gaussian bell formed by the 

284 smearing of the inherent shape of the luminosity peak due to the LumiCal resolution. 

285 The second effect is difficult to precisely correct because of the strong dependence on the 

286 position of the energy cut, and because of the uncertainties of the inherent width of the lumi- 

287 nosity peak and of the energy resolution, as well as the strong correlations between the fitted 
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288 parameters that dominate the spectrum in the peak area (see Eqs. 12 and 13). However, if the 

289 energy cut is made at a sufficient distance from the peak, the second effect becomes negligible, 

290 and the energy-resolution effect can be precisely corrected based on the parametrization of the 

291 functional form of the experimental spectrum after deconvolution of the ISR. 

292 If the inherent width of the luminosity peak is neglected, 3S{Ecm) can be parametrized by the 

293 beta distribution, 

m{EcM) = bo5{EcM-Eo) + I ^i^cm(^o -^cm)"' ;£cm < (^3) 

294 One may recall here that the use of a constant standard deviation a in Eq. 12 is an approxima- 

295 tion, as a depends on the particle energy, and is thus different for different Ecm- The systematic 

296 error induced by the energy resolution of LumiCal can now be expressed as. 



/ ^i^*{EcM) - :nEcM))dEcM 
f ^^iEcM)dEcM 

F , 



(14) 



297 This expression can now be estimated by numerical integration based on the fitted parameters 

298 of .'E*{Ecm) (Eqs. 12 and 13). Even though the reproduction of the integral count by integration 

299 of the fitted function has in principle limited accuracy, rather accurate prediction of the relative 

300 error (Eq. 14) is achieved. The fit was performed on the deconvoluted histogram with the 

301 fixed parameters £0 = 3 TeV and a = 13.7 GeV, while ^o-3 were varying freely. The resulting 

302 parameter values with the associated errors and the correlation matrix are shown below. 



Eq = 3000 
a = 13.7 

Z^o = (676.3 ± 5.3) X 10^ 
Z7i =539±23 
b2 = -0.66±0.15 
^3 = -0.446 ±0.010 

(15) 

0.75 0.64 0.80 \ 

1 0.95 0.98 
0.95 1 0.90 
0.98 0.90 1 / 



1 

0.75 
0.64 
V 0.80 
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(E„-E,J/E^(%) AN/N(10-') 

Figure 6: Left: relative deviation of the peak count induced by the LumiCal energy resolution 
in the reconstructed spectrum as a function of the peak region expressed as fraction 
of the nominal CM energy Eq, compared to the predicted value based on the fitted 
spectrum. Right: histogram of the normalized residual deviations of the peak count 
after correction (see text), for Ecut more than 2.5% away from Eq. 

303 The relative deviation of the count in the reconstructed peak is shown in the left pane of Fig. 

304 6 as a function of the relative distance of the energy cut to the peak energy in percent (black 

305 line). The predicted deviation according to Eq. 14 is also shown for comparison (blue line). 

306 There is an excellent agreement between the predicted and the simulated deviations. To take 

307 a safe distance from the peak, only points for which Ecut is more than 2.5% away from Eq, 

308 corresponding to about 5 a of the fitted peak, will be considered in the following. 

309 The fluctuations of the simulated deviation curve in Fig. 6 are of statistical nature. These 

310 fluctuations can be used as an external measure of the statistical uncertainty of the counting bias 

311 in the simulation. In the right pane in Fig. 6, the histogram of these fluctuations is shown, 

312 calculated as residual deviations after correction, for Ecut more than 2.5% away from Eq. The 

313 RMS of the fluctuations corresponds to a relative statistical uncertainty of 0.24 x 10^^ with 

314 respect to the the peak count in the top 5%. The relative deviation in the top 5% estimated from 

315 Eq. 14 is —0.29 x 10^^. The mean residual bias after correction is (0.05 ±0.03) x 10^^. 

316 Similar procedure was applied to estimate the relative bias and the residual uncertainty in the 

317 tail region from 80% to 90% of Eq. The RMS of the fluctuations is 0.79 x 10"^, the uncorrected 

318 deviation is +0.32 x 10^^^, and the residual deviation after correction is (0.09 it 0.09) x 10^^. 

319 3.5 The Electromagnetic Deflection 

320 To estimate the counting loss due to the EMD, the angular selection was applied once before and 

321 once after the deflection in the simulation, and the relative difference in the resulting number 

322 of events was calculated. The EMD counting loss above 95% of the nominal CM energy is 

323 (-0.50 ± 0.05) X 10"^ . In the tail from 80 to 90% of the nominal CM energy, the EMD counting 

324 loss is (-1.08 ±0.08) X 10"^ 
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325 



4 Conclusions 



326 A method of invariant counting of Bhabha events was presented. The number of Bhabha events 

327 within a given range of scattering angles in the collision frame, and in a given range of Ecm 

328 is reconstructed. The coixesponding limits can be used for the cross-section integration in a 

329 straightforward way. In this way the luminosity expression (Eq. 1) is essentially insensitive to 

330 the beam-beam effects. 

331 The residual systematic uncertainties of the Bhabha count with the presented methods were 

332 estimated by MC simulations. In addition, the systematic uncertainty due to the EMD-induced 

333 counting loss was estimated and found to be small. The residual relative errors in the top 5%, as 

334 well as in the tail from 80 to 90% of the nominal CM energy are hsted in Tab. 2. Beam-beam 

335 effects in the luminosity measurement at CLIC can be corrected with a few permille precision. 

336 The luminosity spectrum from 2.2 TeV to the maximum CM energy can be reconstructed 

337 from the deconvoluted Bhabha-event spectrum (Fig. 5). 

Table 2: Relative residual error after correction of different systematic effects in luminosity mea- 
surement in the peak above 95% and the tail from 80 to 90% of the nominal CM energy. 



# 


Effect 


Top 5% 


80 -90% of £o 






(10-3) 


(10-3) 


1 


Beamstrahlung-ISR angular loss 


-0.1 ±0.4 


-3.6 ±1.8 


2 


High j8„/, ^ 


-0.019±0.008 


-2.7 ±0.1 


3 


ISR energy-loss 


1.3 ±2.0 


-2.3 ±3.9 


4 


Energy resolution 


0.05 ±0.03 


0.09±0.09 


5 


EMD counting loss (uncorrected) 


-0.50 ±0.05 


-1.08±0.08 




Total 


1.4 ±2.0 


4.4 ±4.3 




Total (corrected for #2) ^ 


1.4 ±2.0 


2.7 ±4.3 



The bias under #2 is included in #1. 
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